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Three-dimensional quantum percolation problems are 
studied by analyzing energy level statistics of electrons on 
maximally connected percolating clusters. The quantum per- 
colation threshold p q , which is larger than the classical perco- 
lation threshold p c , becomes smaller when magnetic fields are 
applied, i.e., p q (B = 0) > p q (B 7^ 0) > p c . The critical ex- 
ponents are found to be consistent with the recently obtained 
values of the Anderson model, supporting the conjecture that 
the quantum percolation is classified onto the same univer- 
sality classes of the Anderson transition. Novel critical level 
statistics at the percolation threshold is also reported. 

Studies of percolation [jlj have been attracting much in- 
terest and a number of applications have been proposed 
due to its simplicity. One of its important applications 
is to the transport properties in three-dimensional (3D) 
solids. However, at low temperature, the classical per- 
colation picture is insufficient and the quantum effects 
should be seriously taken into account; thus localization 
due to the quantum interference effect vs. percolation 
becomes a very interesting problem || . 

Even if a maximum cluster percolates through the sys- 
tem, the wave function may be localized and be unable 
to carry current. Thus the quantum percolation thresh- 
old p q is generally larger than the geometrical classical 
percolation threshold p c . At p q , the localization length 
diverges and the wave function extends through the per- 
colating cluster, ft is often debated whether this transi- 
tion is the same as the Anderson transition ||] . If this is 
true, the critical exponent v describing the divergence of 
the localization length in the quantum percolation prob- 
lem would be the same as that in the Anderson transi- 
tion, and the scaling theory of localization would be valid 
[4,5|. Then, a number of studies on Anderson transition 
[ 5, 7| can be used to study quantum percolation problems. 
However, estimates of the critical exponent v often differ 
from that obtained in the Anderson tight binding model. 
For example, renormalization group analyses of the quan- 
tum percolation problem gave v — 2.1 v = 1.9 ±0.5 
@, and v = 1.86 ± 0.05 Q, which are considerably 
larger than the recent estimate v = 1.57 ± 0.02 using 
the Anderson model j[^,[l3), as analyzed by the transfer 
matrix method fjlfl . 

The transfer matrix method, which has been success- 
fully used in the study of the Anderson tight binding 
model, however, cannot be applied to the quantum per- 
colation problem. This is because there is always a fi- 
nite probability that all bonds connecting one plane of 
cross section to the next are vanishing, so that the sys- 
tem breaks into disconnected pieces. To overcome this 



difficulty, it was proposed [jl5| that a missing bond be 
replaced by a bond whose transfer energy is very small 
(~ 10~ 5 ) compared to the transfer energy between the 
connected bonds. A network of percolating wires was 
also considered where the effect of the missing bonds was 
exactly taken into account to calculate the two-terminal 
conductance Jl6| . Though the percolation thresholds de- 
termined in the above studies are consistent with the 
analytical studies of quantum percolation by the Pade 
approximation Jl7| , p^[ , the above numerical study of con- 
ductance and an analytical one |lj| suggest low values of 
v (0.73 and 0.38, respectively) compared to that in the 
Anderson model. 

The analysis of energy levels is free from such geomet- 
rical problems. The Thouless number calculated from 
the energy levels in different boundary conditions gives a 
good estimate of the quantum percolation threshold in a 
three-dimensional system pc|^l[ , but the exponent v is 
1.95 ± 0.12, which is larger than the estimate in the An- 
derson model. Thus, the estimates of critical exponent 
are widely scattered. Even in 2D where all the states are 
believed to be Anderson-localized [Q , some of the above 
studies indicate that there is a localization-delocalization 
transition ||@. 

Another method to extract information on the Ander- 
son transition from eigenvalues is the finite-size scaling 
analysis of energy level spacing |22]-|29]|. Berkovits and 
Avishai applied this method and obtained reliable 
values of p q and v in the quantum percolation problem. 
The obtained value of the exponent v = 1.35 ±0.1 for 
the quantum percolation problem coincides with that 
estimated in the Anderson model by the energy level 
statistics pq , v = 1.45 ±0.1, though it is still slightly 
small compared to that estimated by the transfer matrix 
method for the Anderson tight binding model. 

One of the most striking features of the Anderson tran- 
sition is that it is classified onto a few universality classes 
according to the symmetry under the operation of time 
reversal. Therefore, to study the quantum percolation 
problem in relation to the Anderson transition, the sen- 
sitivity of p q and v to the breaking of the time reversal 
symmetry (TRS) is a very important subject that needs 
to be investigated. 

Here we study the energy level statistics of the quan- 
tum percolation model and discuss the effect of the break- 
ing of TRS. To be specific, we restrict ourselves to the 
bond percolation problems. A peculiar form of the uni- 
versal function of the level spacing distribution at the 
critical point, which is related to the sensitivity to the 
boundary conditions, is also reported. 
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To describe the 3D quantum bond percolation model, 
we consider the following simple Hamiltonian, 



H — yfoijajdj + h.c), 



(1) 



where (ij) denotes the nearest neighbors. The transfer 
integral is defined as 



tij — 



exp(— 2m<pij ;) (for connected bond) , . 

(for disconnected bond) ' 



where the energy unit is the absolute value of the trans- 
fer energy between connected bonds. Bonds are ran- 
domly connected with probabilities p. 4>ij is the Peierls 
phase due to magnetic fields. When all the Peierls phases 
are set to 0, the Hamiltonian is time reversal symmet- 
ric, and we call it the time reversal symmetric (TRS) 
model. When the phases are not vanishing, the Hamil- 
tonian is generally not time reversal symmetric. We set 
— 7r < 4>ij < 7r randomly, and call this situation non- 
TRS model hereafter. The underlying lattice is a three- 
dimensional cube of length L with periodic boundary 
conditions imposed. 

For each realization of connected bonds, we first iden- 
tify the maximally connected percolating cluster, and 
then we diagonalize the Hamiltonian corresponding to 
this cluster by Lanczos method. The calculation is per- 
formed for N different realizations of random bond con- 
figurations, where N = 580,300,175 and 110 for sizes 
L 3 = 12 3 , 15 3 , 18 3 and 21 3 , respectively. These parame- 
ters are chosen so that the number of eigenvalues for each 
system size exceeds 10 6 . 

We analyze the eigenenergies in the region where the 
density of states is smooth (Fig. 1). In the actual simu- 
lation, we take \E\ — 0.2 ~ 0.8. Energy spectra are then 
unfolded |jl 32 1, and the distribution function P(s) of 
adjacent level spacings s is calculated, which is normal- 
ized as 



/>oo />oo 

/ dsP(s) = 1 , / dssP(s) = 1. 
Jq Jo 



(3) 



For a sufficiently large size, we expect the Poissonian 
behavior Pp(s) — exp(— s) for localized states, while it 
should be the Wigner-Dyson type for extended states 
HIH where P W i g (s) °c s /3 exp(-^s 2 ) [fi = 1 in the 
presence of TRS and (3 — 2, otherwise). Varying the 
probability p and the size L, we extract the information 
of the quantum percolation transition. 

Plots of P(s) as a function of the bond occupation 
probability for the TRS model and the non-TRS model 
(L = 21) are displayed in Fig. |2| It is clearly seen 
that the expected crossover from a Poissonian behavior 
to that of Wigner-Dyson type is manifested as we increase 
the value of p. 

In order to obtain the critical value of the probability 
p q and the critical exponent for the divergence of the 
localization length v, we define I(s) and A(p, L) as p4j 



and 



A(p,i) 



/(*)= [ S p(s')ds', 
Jo 



Jo" tjsjds - f°° I P (s)ds 
j;°I Wlg ( S )d S -j; o I P ( S )ds- 



(4) 



(5) 



Here we set so = 1-2. We also use other quantities 
to characterize the critical behavior of level statistics, 
PH and have confirmed that the results of the transition 
points and the critical exponents are almost the same. 

Denoting £(p) as the localization length or the corre- 
lation length which diverges as 
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b-p q r 

this function is expected to behave as 

A(p,L) = /[L/£(p)] 

= a Q + a 1 [L{p-p q y\ 1 l v 

a 2 (lL(p-p q y}^) 2 + 



(6) 



(7) 



near the critical probability p q . 

In Fig. |[ we show A(p, L) as a function of p for different 
sample sizes L = 12, 15, 18 and 21. All curves cross at a 
single point. By fitting A(p, L) to eq.(0), we can estimate 
the critical exponent as well as p q . 

In Table I, we summarize our results of the \ 2 fit 
to eq.([7|). Here, the error bars of raw data as well as 
the fitting parameters are estimated using the bootstrap 
method (l^ |l^,|34||35j . Sizes less than L — 9 give an un- 
reasonably large ;p7 so we have excluded them. 

The earlier result by Berkovits and Avishai |3(| for the 
TRS model gives 1.35 ± 0.1, which is smaller than the 
transfer matrix result for the Anderson model, 1.57±0.02 
J36|. Here we have considered two possibilities: (i) the 
effect of energy levels of small disconected clusters, and 
(ii) the corrections to scaling due to small sizes. 

Let us begin with the first possibility. Including the 
levels of small disconnected clusters |p0|| , in principle, in- 
fluences the results, since levels in the isolated clusters 
are not correlated with those in the percolating cluster, 
which destroys the level correlation. However, the results 
from all the energy levels (TRS (all)) almost coincide with 
those extracted from the information of the eigenenergies 
only for the states on the percolating cluster (TRS) . This 
might be due to the fact that we are observing the region 
far from the classical percolation threshold p c « 0.249, 
and most of the states belong to the percolating clusters. 
|5?J We conclude that the possibility (i) is irrelevant, and 
thus discuss next possibility (ii). 

The analysis by finite-size scaling of the level statistics 
tends to give smaller values [B2 24 2^|3(J compared to 
that by the transfer matrix method combined with finite- 
size scaling R,[l2|. It is possible that such discrepancies 
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arise from the small size of the system, where corrections 
to scaling are not negligible J38|— pp| , p~3| | . To investigate 
the corrections to scaling, we omit the data for smaller 
system sizes In the TRS model, excluding smaller 
system sizes tend to increase v, which saturates around 
v w 1.5, in agreement with the Anderson model. 

Breaking the time reversal symmetry (non-TRS 
model) decreases the value of the transition point p q Q . 
This can be attributed to the destruction of the con- 
structive quantum interference among the time-reversed 
paths, which supports the idea that this transition is in- 
deed the Anderson transition. In the non-TRS model, 
v tends to increase by excluding smaller system sizes, 
but the fitting becomes unstable when only the data for 
L = 18 and 21 are used. We conclude that the value of v 
in the non-TRS model is consistent with that of the An- 
derson mdocl belonging to the unitary universality class 
Hi (i/= 1.43 ±0.02). 

Now we focus our attention on the level statistics at 
the critical point p q . In Fig. |[ we plot the distribution 
functions P(s) for the nearest-neighboring level spacing 
s at Pq for various system sizes. They seem to be in- 
dependent of the sizes, but are different from that for 
the Anderson model with the periodic boundary condi- 
tion (p.b.c.) and from that with the fixed boundary con- 
dition (f.b.c). Recently, it was reported that the form 
of -P(s) at the critical point changes with the change of 
boundary conditions p^ j. The critical distribution func- 
tion P(s) obtained in the present study lies between the 
two extreme cases. This is consistent with the fact that 
some of the bonds are missing at the edges of the system 
in the percolation problem, so even if we adopt the peri- 
odic boundary conditions, the model is effectively closer 
to the one with f.b.c. Here, we have plotted the results 
for the non-TRS model, but the same behavior is also 
observed in the TRS model. 

In conclusion, we have analyzed the three-dimensional 
percolation problems by the finite-size scaling analysis 
for the energy levels. Breaking the time reversal sym- 
metry lowers the quantum percolation threshold. Both 
the time reversal symmetric and non-symmetric models 
exhibit the critical exponent v consistent with the results 
of Anderson tight binding model. At the critical point, 
the energy level spacing distribution, which characterizes 
the level statistics, becomes size-independent, which dif- 
fers from the one in the Anderson model with p.b.c. or 
f.b.c. 

At present, we cannot clearly show that the TRS model 
and the non-TRS model belong to different universality 
classes. It took large-scale numerical effort to distinguish 
one universality class from the other even for the Ander- 
son model Jl^j , since highly accurate data of less than 
0.2% relative error are required. This high accuracy is 
very difficult to achieve by level statistics for the quan- 
tum percolation problem. Much more numerical effort is 
necessary for this purpose. 
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FIG. 1. Normalized density of states of the quantum per- 
colation problem for the non-TRS model. The solid line is 
for the maximally connected cluster, while the dashed line 
corresponds to all the energy levels including small clusters. 

FIG. 2. Level spacing distribution P(s) for the TRS model 
(a) and the non-TRS model(b) for L = 21. With an increase 
in p, the transition from a Poisson distribution to that of 
Gaussian orthogonal ensemble is seen in Fig. 2(a), while the 
transition from a Poisson distribution to that of Gaussian 
unitary ensemble is seen in Fig. 2(b). 

FIG. 3. A(p, L) as a function of p for different sample sizes 
for the TRS model (a) and the non-TRS model(b) obtained 
from the energy of the percolating cluster. Energy range is 
\E\ = 0.2 ~ 0.8. At p q , all curves cross at a point where 
A(p,L) becomes scale-independent. Bars around the data 
points indicate error bars estimated by the bootstrap method, 
and the curves are the results of the fit to eq.(Q). Insets are 
the scaling plots. 

FIG. 4. Critical level spacing distribution for the non-TRS 
model, which shows intermediate behavior between that of 
Anderson model with the p.b.c. (solid line) and those with 
f.b.c. (dashed line). 

TABLE I. The % 2 fitting of critical point p q and critical 
exponent v with their standard deviations. The values of the 
TRS model and the non-TRS model are only for percolating 
cluster, while TRS (all) is obtained from all energy levels. 
TRS (all,BA) is the result reported by reported Berkovits 
and Avishai [30]. 
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Pq 


V 




TRS (all,BA) 


L = 


7,9,11,13,15 


0.33 ± .01 


1.35 ± 


.10 


TRS (all) 


L = 


12,15,18,21 


0.321 ± .001 


1.48 ± 


.10 


TRS 


L = 


12,15,18,21 


0.324 ± .001 


1.46 ± 
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TRS 
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15,18,21 


0.323 ± .001 
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1.50 ± 
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non-TRS 


L = 


12,15,18,21 


0.309 ± .001 


1.24 ± 


.08 


non-TRS 


L = 


15,18,21 


0.308 ± .001 


1.44 ± 


.16 
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